SPINODAL INSTABILITIES IN NUCLEAR MATTER IN A 
STOCHASTIC RELATIVISTIC MEAN-FIELD APPROACH 

S. Ayik 1 ]] 0. Yilmaz 2 , N. Er 2 , A. Gokalp 2 , and P. Ring 3 
1 Physics Department, Tennessee Technological University, Cookeville, TN 38505, USA 
2 Physics Department, Middle East Technical University, 06531 Ankara, Turkey 
3 Physics Department, TU Munich, D-851^8 Garching, Germany 

(Dated: July 10, 2009) 

Abstract 

Spinodal instabilities and early growth of baryon density fluctuations in symmetric nuclear mat- 
ter are investigated in the basis of stochastic extension of relativistic mean-field approach in the 
semi-classical approximation. Calculations are compared with the results of non-relativistic calcu- 
lations based on Skyrme-type effective interactions under similar conditions. A qualitative differ- 
ence appears in the unstable response of the system: the system exhibits most unstable behavior 
at higher baryon densities around pb = 0.4 pq in the relativistic approach while most unstable 
behavior occurs at lower baryon densities around pb = 0.2 po in the non-relativistic calculations 
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I. INTRODUCTION 



Spinodal instability provides a possible dynamical mechanism for fragmentation of a hot 
piece of nuclear matter produced in heavy-ion collisions. Small amplitude density fluctua- 
tions grow rapidly and lead to break-up of the system into an ensemble of clusters In 
coming years, experimental investigations of multi-fragmentation reactions in neutron rich 
nuclear system will provide further understanding of isospin dependence of nuclear matter 
equation of state at low densities. In theoretical side, extensive investigations of spinodal in- 



stabilities have been carried out in the basis of stochastic transport models 



, y, Q, Q, E 



6|. In 



particular, the recently proposed stochastic mean-field approach provides a useful tool for a 
description of dynamics of density fluctuations in the spinodal region 7|. It has been demon- 
strated that the stochastic mean-field approach incorporates the one-body dissipation and 
the associated fluctuation mechanism in accordance with the quantal-dissipation fluctuation 
relation. The approach gives rise to the same result for dispersion of one-body observables 
that was obtained in a variational approach in a previous work 8j. Furthermore, in recent 
studies (9), n]] by projecting onto macroscopic variables, we deduce transport coefficients 
for energy dissipation and nucleon exchange in low-energy heavy-ion collisions, which have 
;he similar form with those familiar from the phenomenological nucleon exchange model 



These investigations provide a strong support for the fact that the stochastic mean- 



field approach is a powerful tool for describing low energy nuclear collisions and spinodal 
dynamics. 

In a recent work, we studied the early development of spinodal dynamics of nuclear matter 
in the basis of the stochastic mean-field approach by employing density- dependent Skyrme- 
type effective interactions 121 ] . In the present work, we carry out a similar investigation of 



early development of density fluctuations in spinodal region of nuclear matter by employing 
the stochastic extension of the relativistic mean-field theory ljj, [l^] . It has been shown in 
recent years that the nuclear many-body system is in principal a relativistic system driven 
by dynamics of large relativistic attractive scalar and repulsive vector fields. Both fields 
are not much smaller than the nucleon mass and therefore the average nuclear field should 
be described by Dirac equation. For large components of Dirac spinors, two fields nearly 
cancel each other leading to relatively small attractive mean field. The small components 
add up leading to a very large spin orbit term, which is known since early days of nuclear 
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physics. Relativistic models have been used with great success to describe nuclear structure. 
In recent years, the approach has also been applied for description of nuclear d yna mics 



extended in the framework of time- dependent covariant density functional theory 15|, [16]. 



A number of investigations have been carried out on spinodal instabilities in nuclear matter 



employing relativistic mean- field approaches 17|, [l8|, [l9|. In this work, we consider the 



stochastic extension of the relativistic mean-field theory in the semi-classical approximation. 
As illustrated in the non- relativistic limit, stochastic extension of the mean- field theory 
provides a powerful approach for investigating dynamics of density fluctuations. Employing 
the stochastic extension of the relativistic mean-field approach, we investigate not only 
spinodal instabilities but also the early development of density fluctuations in symmetric 
nuclear matter. 

In Section 2, we briefly describe the stochastic extension of the relativistic mean- field 
theory in the semi-classical approximation. In Section 3, we calculate early growth of baryon 
density fluctuations, growth rates and phase diagram of dominant modes in symmetric 
nuclear matter. Conclusions are given in Section 4. 

II. STOCHASTIC RELATIVISTIC MEAN-FIELD THEORY 



The stochastic mean-field approach is based on a very appealing stochastic mode 



21 



pro- 
22) 



posed for describing deep-inelastic heavy-ion collisions and sub-barrier fusion [20 
In that model, dynamics of relative motion is coupled to collective surface modes of collid- 
ing ions and treated in a classical framework. The initial quantum zero point and thermal 
fluctuations are incorporated into the calculations in a stochastic manner by generating an 
ensemble of events according to the initial distribution of collective modes. In the mean-field 
evolution, coupling of relative motion with all other collective and non-collective modes are 
automatically taken into account. In the stochastic extension of the mean-field approach, the 
zero point and thermal fluctuations of the initial state are taken into account in a stochastic 
manner, in a similar manner presented in refs. {20, 21, 2^]. The initial fluctuations, which 
are specified by a specific Gaussian random ensemble, are simulated by considering evolu- 
tion of an ensemble of single-particle density matrices. It is possible to incorporate quantal 
and thermal fluctuations of the initial state into the relativistic mean-field description in a 
similar manner. 
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In refs. |23j, |24|, the authors derived a relativistic Vlasov equation from the Walecka 
model in the local density and the semi-classical approximation. In the Walecka model, 
interaction between nucleons are mediated by a scalar meson with mass m s and a vector 
meson with mass my, with respective fields denoted as <f> and V^. Introducing phase space 
distribution function f(f,p,t) for the nucleons, following relativistic Vlasov equation has 
been obtained , 

d 



Ot 



f{f,p,t) +v- V r f(f,p,t) - V r h(r,p) ■ V p /(r,p,t) = 0. 



(1) 



where v = p*/e* and h = e* + gyVo. The coupling constants of the mesons and the nucleon 
are denoted by g s and gy, for the scalar and the vector mesons, respectively. In these 
expressions, p* = p — g v V and e* = (p* 2 + M* 2 ) 1 / 2 with M* = M — g S (j). The nucleon mass 
is denoted by M. In the mean-field approximation, the meson fields are treated as classical 
fields and their evolutions are determined by the field equations, 



dt 2 



^ - V 2 + m 2 



(r,t) = gsPs{r,t) 



(2) 



and 



dt 2 



- V 2 + m 2 



V v (r,t) = g vPv {r,t). (3) 

In these expressions, the baryon density p (r,t) = pu(r,t), the scalar density p s (f,t), and 
the current density py(f,t) can be expressed in terms of phase-space distribution function 
as follows, 

d 3 p 



Pb^) = 7 



and 



Ps(r,t) = 7 



py{r,t) = 7 



(27T) 

d 3 p M* 
(2tt) 3 ~ 



rf(r,P,t), 



f(r,P,t), 



d 3 p p 



(27T) 



-f(r,P,t) 



(4) 
(5) 

(6) 



where 7 = 4 is the spin-isospin degeneracy factor. The original Walecka model gives a 
nuclear compressibility that is much larger than the one extracted from the giant monopole 
resonances in nuclei. It also leads to an effective nucleon mass which is smaller than the 
value determined from the analysis of nucleon-nucleon scattering. In order to have a model 
which allows different values of nuclear compressibility and the nucleon effective mass, it is 
possible to improve the Walecka model by including the self-interaction of the scalar mesons 



or by considering density dependent coupling constants. However, in the present exploratory 
work, we employ the original Walecka model without including the self interaction of the 
scalar meson. 

In the stochastic mean-field approach an ensemble {f x {r,p,t)} of the phase-space distri- 
butions is generated in accordance with the initial fluctuations, where A indicates the event 
label. In the following for simplicity of notation, since equations of motions do not change 
in the stochastic evolution, we do not use the event label A for the phase-space distributions 
and also on the other quantities. However it is understood that the phase-space distribution, 
scalar meson and vector meson fields are fluctuating quantities. Each member of the ensem- 
ble of phase-space distributions evolves by the same Vlasov [l] equation according to its own 
self-consistent mean-field, but with different initial conditions. The main assumption of the 
approach in the semi- classical representation is the following: In each phase-space cell, the 
initial phase-space distribution f(f,p, 0) is a Gaussian random number with its mean value 



determined by f(r,p,0) = fo(r,p), and its second moment is determined by 



f(r,p,0)f(r*,i?,0) = (27r) 6 8(r- f*)5(p- $)f o (r,0)[l - f o (r,0)] (7) 

where the overline represents the ensemble averaging and /o (r , p) denotes the average phase- 
space distribution describing the initial state. In the special case of a homogenous initial 
state, it is given by the Fermi-Dirac distribution fo(p) = l/[exp(eg — Po)/T + 1]. In this 
expression pi = fi — (gy / my)" 2 where po is the chemical potential and p^ is the baryon 
density in the homogenous initial state. 

In this work, we investigate the early growth of density fluctuations in the spinodal region 
in symmetric nuclear matter. For this purpose, it is sufficient to consider the linear response 
treatment of dynamical evolution. The small amplitude fluctuations of the phase-space 
distribution 5f(f,p,t) = f(r,p,t) — fo(p) around an equilibrium state fo(p) are determined 
by the linearized Vlasov equation, 

^6f(r,p,t)+v ■ V r 6f(r,p,t) - V r 5h(r,p,t) • V„/ (p) = 0. (8) 



In these expression the local velocity is v = pje% with = ^p 1 + Mq 2 , Mq = M — gs<Po, 
and small fluctuations of mean-field Hamiltonian is given by, 

5h(f,p,t) = -MjLgrftftQ+gySVoftt) - 5V(f,t) (9) 

e o e o 



The small fluctuations of the scalar and vector mesons are determined by the linearized field 
equations, 

64>(r,t) = g 8 5p s {f,t) (10) 



Q2 

|^-V 2 + m§ 

and 



r%2 

— - V 2 + m\ 5V v (r, t) = gySfrfr *)■ < 11 > 



III. EARLY GROWTH OF DENSITY FLUCTUATIONS 
A. Spinodal Instabilities 

In this section, we employ the stochastic relativistic mean-field approach in small ampli- 
tude limit to investigate spinodal instabilities in symmetric nuclear matter. We can obtain 
the solution of linear response equations (ITl)- (llll) by employing the standard method of one- 
sided Fourier transform in time [25|]. It is also convenient to introduce the Fourier transform 
of the phase-space distribution in space, 

~ —> foe roo — 

6f(k,p,u) = dte™ d 3 re~ ik - r f(r,p,t). (12) 

JO J-oo 

This leads to, 

of{k,p,uj) = —g s 5p s {k,uj) - g v 5p h (k,uj) + g v — ■ dp y (k,u) ^ — f + i — f 

\ e o e o / to — v ■ k uj — Vo • k 

(13) 

where 5f(k,p,0) denotes the Fourier transform of the initial fluctuations, and we use the 
short hand notation, = g1/{k 2 + mg), g^ = g\j(k 2 + rn\). In this expression, the 
fluctuations of the meson fields are expressed in terms of Fourier transforms of the scalar 
density Sp s (r, t), the baryon density Spb(r, t) and the current density 5pv(r, t) fluctuations by 
employing the field equations (pTJ|) - (TTTj) . In Eq. ([13]) only the initial fluctuations of the phase- 
space distribution 5f(k,p,0) is kept, but the initial fluctuations associated with the scalar 
and the vector fields are neglected. In the spinodal region since it is expected to have a small 
contribution, we neglect the frequency terms in the propagators, i.e., — uo 2 + k 2 +m§ ~ k 2 +rrig 
and — uj 2 + k 2 + m\ ~ k 2 + my. Small fluctuations of the baryon density, the scalar density 
and the current density are related to the fluctuation of phase-space distribution function 
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5f(k,p,uj) according to, 
5p h (k,u) = 7 ' 

5p s (k,u) 



7 



7 



(2tt) 3 
d 3 p 
(2^)3 
d 3 p 



Sf(k,p,uj 
5' M 



(14) 



M 



/o(p) + 



M„* 



^v9fP'^v(^^) ~ g&$Ps(k,u) ) foip) + A -^-8f(k,p,uj) 



(15) 



and 



5p v (k,u) 
d 3 p 



7 



7 



(27T) 



d 3 p 
(2^)3 

P 



P 



/o(p) + ^-5f{k,p,uj) 



#v— P • <Vv( fc > w ) - 2v ; 

e e 



^s-§-P s Ps(k,uj) I /o(p) + ^Sf(k,p,uj) 

e / e 



p 



(16) 



(17) 



Multiplying both sides of Eq. (fT3]) by Mq/cq, 1, p/eg and integrating over the momentum, 
we deduce a set of coupled algebraic equations for the small fluctuations of the scalar density, 
the baryon density and the current density, which can be put in to a matrix form. Here we 
investigate spinodal dynamics of the longitudinal unstable modes. For longitudinal modes 
the current density oscillates along the direction of propagation, 5p v (k,tu) = 8py(k,u)k. 
Then, for the longitudinal modes, the set of equations become, 

< A x A 2 A 3 \( 8p v (k,oo)\ 

B 1 B 2 B 3 5p s (k,iu) 
\ Ci C 2 C 3 J \5p h (k,uj) J 
where the element of the coefficient matrix are defined according to, 

B x B 2 B 3 = -fjyXv&w) l + ^|xs(fc,w) +g v Xs{k,uj) 
\C Y C 2 C 3 ) \l + iftxh(J*> w ) ~9ixv{k,uj) +g$xv(k,u) J 
In this expression, Xb(^> Co ')> Xs(k, ui) and xv(k,u) denote the long wavelength limit of rela- 
tivistic Lindhard functions associated with baryon, scalar and current density distribution 
functions, 



f S h (k,u)\ 
S s (k,uj) 

\ Sy(k,Uj) J 



(18) 



Xs(k,u) 



7 



d 3 p 



p-k/e*\ 

i 



fc • V p /o(p) 
uj — v • k 



(19) 



and the stochastic source terms are determined by 



S s (k,uj) 

\ Sy(k,u) J 



( 



7 



d 3 p 



\ 



M*/el 
\p-k/el ) 



Sf(k,p,Q) 
oj — v ■ k 



(20) 



Other three elements of the coefficient matrix in Eq. fll8p are given by, 



Xs{k,u) = 7 



d 3 p 
(2^)3 



P 2 f ,s MS 2 k ■ V p f (p) 
eo e u — v • k 



Xv(k,uj) = 7 / -^S^P- k 



and 



-» f d 3 p 



(2tt) 3 
ef - (p- kf 



M*k-V p f (p) 



e o 2 to — v • k 



P *3 



hi?) 



{p-k) 2 k-V p h{p) 



e 



(21) 



(22) 



(23) 



u — vo • k 

We obtain the solutions by inverting the algebraic matrix equation, which gives for the 
baryon density fluctuations, 



5p B (k,u) = i 



D 1 S h (k,u) + D 2 S s (k,m) + D 3 S v (k,cu) 
e(k,u) 



(24) 



where D\ = B\C 2 — B 2 C±, D 2 = C\A 2 — C 2 A\ and D 3 = A\B 2 — A 2 B\ and the quantity 
e(k,u) = A3D1 + B 3 D 2 + C3D3 denotes the susceptibility. 

The evolution in time is determined by taking the inverse Fourier transformation in time, 
which can be calculated with the help of residue theorem 24|. Keeping only the growing 
and decaying collective poles, we find, 



6p h (k,t) = 6pUk)e +r * t + 6 P 7(k) 



(25) 



Here, the amplitudes of baryon density fluctuations associated with the growing and decaying 
modes at the initial instant are given by, 



6pl(k) = - 



DiS h (k, u) + £ 2 £s(fc, uj) + D 3 Sv(k, u) 
de(k, u)/duj 



(26) 



The growth and decay rates of the modes are obtained from the dispersion relations, 
e(k, u) = 0, i.e. from the roots of susceptibility. Solutions for the scalar density fluctu- 
ations 5ps{k,uj) and the current density 5py(k,u>) fluctuations can be expressed in a similar 
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FIG. 1: Growth rates of unstable modes as a function of wave numbers in the spinodal region 
at baryon densities p^ = 0.2 po and p^ = 0.4 po at temperature T = 5MeV, (a) relativistic 
calculations, and (b) non-relativistic calculations. 

manner. In the original Walecka model, there are four free parameters, coupling constants 
and meson masses. The binding energy per nucleon at saturation density determines the ra- 
tios of coupling constants to masses. The standard values of the ratios gy(M/my) 2 = 273.8 
and _q| (M/m,s) 2 = 357.4 give binding energy per nucleon 15.75 MeV at saturation density 



m 



14j . These ratios lead to an effective nucleon mass Mq = 0.541M and a compressibility 
of 540 MeV at the saturation density. In numerical calculations, we take for the vector 
meson mass my = 783 MeV, and for the scalar meson mass, mg = 500 MeV. As an ex- 
ample, the upper panel in Fig. [T] shows the growth rates of unstable modes as a function of 
wave number in the spinodal region corresponding to the initial baryon density p^ = 0.2 po 
and p^ = 0.4 p at a temperature T = 5 MeV. The lower panel of Fig. [1] illustrates 
the dispersion relations obtained in the non-relativistic approach with an effective Skyrme 
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FIG. 2: Growth rates of the most unstable modes as function of baryon density in spinodal region at 
temperature T = 5MeV in relativistic calculations (solid line) and in non-relativistic calculations 
(dashed line). 



force 12]. Although direct comparison of these calculations is rather difficult, we observe 
there are qualitative differences in both calculations. The range of most unstable modes 
in relativistic calculations is concentrated around k = 0.6 frrT 1 in both densities, while 
most unstable modes shift towards larger wave numbers around k = 0.8 fm^ 1 at density 
= 0.2 po towards smaller wave numbers around k = 0.5 fmT 1 at density p^ = 0.4 po. 
Growth rates of most unstable modes at density p^ = 0.4 po in relativistic calculations are 
nearly factor of two larger than those results obtained in the non-relativistic calculations, 
while at low density p^ = 0.2 p the growth rates are smaller in relativistic calculations. 
Fig. [2] illustrates growth rates of the most unstable modes as a function of density in both 
relativistic and non-relativistic approaches. We observe the qualitative difference in the un- 
stable response of the system: the system exhibits most unstable behavior at higher densities 
around p^ = 0.4 po in the relativistic approach while most unstable behavior occurs in the 
non-relativistic calculations at lower densities around p^ = 0.2 po- As an example of phase 
diagrams, Fig. [3] shows the boundary of spinodal region for the unstable mode of wavelength 
A = 9.0 fm. Again, we observe that the unstable behavior shifts towards higher densities 
in relativistic calculations. 
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FIG. 3: Boundary of spinodal region in baryon density-temperature plane for the unstable mode 
with wavelengths A = 9 fm in relativistic calculations (solid line) and in non-relativistic calculations 
(dashed line). 

B. Growth of Density fluctuations 

In this section, we calculate the early growth of baryon density fluctuations in nuclear 
matter. Spectral intensity of density correlation function (A;, t) is related to the variance 
of Fourier transform of baryon density fluctuation according to, 



a hh (k,t)(2n) 3 6(k - k') = 5p h (k,t)5~p* h (k',t) (27) 

We calculate the spectral function using the solution fl25|) and the expression (jTj) for the 
initial fluctuations to give, 

a hh (k,t) = (e 2r ^ + e" 2 ^) + b ^ '- (28) 

where 



E*{k) = \D 1 \' 2 K^ 1 + \D 2 \ 2 K^t\D 3 \ 2 K^ + 2D 1 D 2 K^ (29) 



with 



^i^^^f.m-mi (30, 

11 



K , 2 = 7 2 / ^ (<] (?.-*)' f m _ fM (31) 



.2 

2 y d 3 p /p-fc\ r|^(v -A;) 2 



= ' / iVJ Fi + ^*-)V /a(p)[1 " /o(p)1, (32) 

and 

*S = 7 2 / F -j^S / (p)[l - / (p)]. (33) 

Upper and lower panels of Fig. H]sliow the spectral intensity of the baryon density correlation 
function as a function of wave number at times t = 0, t = 20 fm/c and t = 40 fm/c at 
temperature T = 5 MeV in relativistic calculations at densities = 0.2 p an d = 0.4 p 0) 
respectively. We observe that the largest growth occurs over the range of wave numbers 
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FIG. 4: Spectral intensity a^(k,t) of baryon density correlation function as a function of wave 
number at times t = 0, t = 20 fm/c and t = 40 fm/c at temperature T = 5 MeV in relativistic 
calculations at density (a) = 0.2 po and (b) p^, = 0.4 p^. 
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FIG. 5: Same as figured] in non-relativistic calculations. 

corresponding to the range of dominant unstable modes. Spectral intensity in the vicinity of 
most unstable modes of k — 0.6 fm" 1 grows about a factor of ten at density p^ = 0.2 p and 
about a factor of six at density p^ = 0.4 po during the time interval of t = 40 fm/c. Fig. [5] 
shows the similar information calculated in non-relativistic approaches. We notice that at 
density p^ = 0.2 p the behavior of spectral intensity is rather similar in relativistic and non- 
relativistic approache. However, at higher density p^ = 0.4 p , the spectral intensity grows 
slower in the non-relativistic calculations than those obtained in the relativistic approach. 
We note that in determining time evolution Sp^(k,t) with the help of the residue theorem, 
there are other contributions arising form the non-collective pole of the susceptibility e(k,u) 
and from the poles of source terms Sy(k,u), Ss(k,u>) and S-^(k,uj). These contributions, in 
particular towards the short wavelengths, i.e. towards higher wave numbers, are important 
at the initial stage, however they damp out in a short time interval |25] . Since, we do not 
include effects from non-collective poles, we terminate the spectral in Fig. [5] at a cut-off 
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wave number k c rs 0.7 fm 1 — 0.8 fm 1 . Consequently, the expression (128]) provides a 
good approximation for dbb{k,t) in the long wavelength regime below k c . 

Local baryon density fluctuations 5p b (r, t) are determined by the Fourier transform of 
5p^(k,t). Equal time correlation function of baryon density fluctuations as a function of 
distance between two space locations can be expressed in terms of the spectral intensity as, 

a hh (\r-r>\,t) = 5nb(f, t)<yn b (F, t) = J -0 e ^> bb (^t). (34) 

The baryon density correlation function carries useful information about the unstable dy- 
namics of the matter in the spinodal region. As an example, the upper and lower panels of 
Fig. E] illustrates the baryon density correlation function as a function distance between two 
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FIG. 6: Baryon density correlation function a^(x,t) as a function of distance x = \r — r"| between 
two space points at times t = 0, t = 20 fm/c and t = 40 fm/c at temperature T = 5 MeV in 
relativistic calculations at density (a) = 0.2 po and (b) p-^ = 0.4 po- 
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space points at times t = 0, t = 20 fm/c, and t = 40 fm/c at temperature T = 5 MeV in 
relativistic calculations at densities p^ = 0.4 p and p^ = 0.2 p , respectively. Complemen- 
tary to the dispersion relation, correlation length of baryon density fluctuations provides an 
additional measure for the size of the primary fragmentation pattern. We can estimate the 
correlations length of baryon density fluctuations as the width of the correlation function at 
at half maximum. From the figure, we estimate that the correlation length is about the same 
at both densities and temperatures around 3.0 fm, which is consistent with the dispersion 
relation presented in Fig. [TJ Baryon density fluctuations grow faster at p^ = 0.4 p than 
P k = 0.2 pq. Fig. [7] shows the similar information calculated in the non-relativistic approach 



121 ] . The correlation length is around 3.0 fm at p^ = 0.4 p and 3.0 fm at the lower density 
p^ = 0.2 pq. However, unlike the relativistic calculations, the baryon density fluctuations 
grow faster at lower density p^ = 0.2 p than at p^ = 0.4 p , which is consistent result 
presented in Fig. [2] 
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FIG. 7: Same as figure [6] in non-relativistic calculations. 
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IV. CONCLUSIONS 



12 1 that the stochastic mean- field 



It has been demonstrated in recent publications [7l.l9l.flO. 
approach incorporates both the one-body dissipation and the associated fluctuation mech- 
anism in a manner consistent with the fluctuation-dissipation theorem of non-equilibrium 
statistical mechanics. Therefore the approach provides a powerful tool for investigating 
dynamics of density fluctuations in low-energy nuclear collisions. In a similar manner, it 
is possible to develop an extension of the relativistic mean-field theory by incorporating 
the initial quantal zero point fluctuations and thermal fluctuations of density in a stochastic 
manner. In this work, by employing the stochastic extension of the relativistic mean-field ap- 
proach, we investigate spinodal instabilities in symmetric nuclear matter in the semi-classical 
framework. We determine the growth rates of unstable collective modes at different initial 
densities and temperatures. Stochastic approach also allows us to calculate early develop- 
ment of baryon density correlation functions in spinodal region, which provides valuable 
complementary information about the emerging fragmentation pattern of the system. We 
compare the results with those obtained in non-relativistic calculations under similar condi- 
tions. Our calculations indicate a qualitative difference in behavior in the unstable response 
of the system. In the relativistic approach, the system exhibits most unstable behavior at 
higher baryon densities around = 0.4 p , while in the non-relativistic calculations most 
unstable behavior occurs at lower baryon densities around = 0.2 p . In the present 
exploratory work, we employ the original Walacka model without self-interaction of scalar 
meson. The qualitative difference in the unstable behavior may be partly due to the fact 
that the original Walecka model leads to a relatively small value of nucleon effective mass 
of M* = 0.541M and a large nuclear compressibility of 540 MeV. On the other hand, the 
Skyrme interaction that we employ in non-relativistic calculations gives rise to a compress- 
ibility of 201 MeV [121 ] . It will be interesting to carry out further investigations of spinodal 
dynamics in symmetric and charge asymmetric nuclear matter by including self-interaction 
of the scalar meson and also including the rho meson in the calculations. Inclusion of 
the self-interaction of scalar meson allows us to investigate spinodal dynamics over a wide 
range of nuclear compressibility and nuclear effective mass. We also note by working in 
the semi-classical framework, we neglect the quantum statistical effects on the baryon den- 
sity correlation function, which become important at lower temperatures and also at lower 
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densities. 
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